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WTO -DIMENSIONAL COMPRESSIBLE BLOW IN CONICAL 
MIXED -FLCW COMPRESSORS 
By John D. Stanitz 


SUMMARY 

A general method of analysis is developed for the two- 
dimensional, steady, compressible flow through mixed -flow compres- 
sors in which the center line of the passage generates a right cir- 
cular cone about the axis of the compressor. Tie two-dimensional, 
radial -discharge compressor is a special case in which the cone 
angle is equal to 180°. Tie variables taken into account Eire: 

(1) impeller tip Mach number, (2) compressor flow rate, (3) blade 
shape (curvature), (4) passage height, (5) number of blades, and 
(6) compressor cone angle. Relaxation methods are used to solve 
the resulting nonlinear differential equation for the stream 
function. Special attention is paid to logarithmic-spiral blades, 
of which the straight blade (lying on a conic radius) 1 b a partic- 
ular case. 

As a result of the analysis, it iB concluded that the solution 
obtained for a given cone angle also applies to certain other cone 
angles (that is, other mixed-flow compressors) with a fewer or a 
greater number of passages but with the same included passage angle, 
and so forth. It is also concluded that mixed-flow compressors 
with the same number of flow passages as radial-discharge com- 
pressors and therefore with smaller included passage angles have 
lower peak blade loadings and lower maximum relative velocities 
than the corresponding radial -discharge compressors. 


The general analysis also applies to inward -flow turbines. 

In fact, the solution obtained for a centrifugal compressor with 
smooth (shockless) entry is also the solution (with the flow direc- 
tion and rotation reversed) for an inward-flow turbine with the 
same design characteristics (that is, the same rotor) and with 
shockless entry. 
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A numerical example is presented consisting of a radial- 
discharge compressor (cone angle equal to l 80 °) with constant flow 
area and 20 straight, radial blades operating at a tip Mach number 
of 1.5. The results of this example are given by plots of the 
streamlines, lines of constant pressure ratio, and lines of constant 
Mach number. For the conditions of this example, a wheel-type eddy 
forms on the driving face of the blade; the velocities .and pres- 
sures at the impeller tip are reasonably uniform; any nonuniformity 
in the flow leaving the impeller tip adjusts itself rapidly; the 
maximum local Mach number is 0.64 and occurs along the trailing 
face of the blade at 68 percent of the tip radius; and the computed 
value of the slip factor is 0.90. 

* 

In addition, a simplified analysis for straight blades lying 
on conic radii is presented that can be used to determine the 
streamlines, the pressure distribution, and the velocity profiles 
within the Impeller except near the tip (and inlet). 

IHTB0DDCTI0IT 

Increased knowledge of flow conditions within centrifugal com- 
pressors can result in improved compressor performance . For 
example, boundary- layer growth and flow separation, which affect 
the compressor efficiency, can be controlled by suitable changes 
in the compressor design provided the effects of these changes on 
velocity and pressure gradients within the compressor are known. 

For a given set of operating conditions, the flow conditions 
within centrifugal compressors depend upon the geometry of the com- 
pressor (three-dimensional effects) and upon the properties of the 
fluid (compressibility and viscosity). Most treatments of the 
problem up to the present time have been concerned with the two- 
dimensional-flaw effects f oir incompressible, nonviscous fluids (for 
example, references 1 to 5). 

In the analysis reported herein, compressibility is considered, 
which is especially important in centrifugal compressors because 
the large pressure ratios per stage result in density changes that 
affect the fluid velocities and therefore the streamlines, the 
pressure gradients, and so forth. In addition, compressibility is 
important in regions of supersonic flow where shock phenomena may 
develop. A method is developed for determining the two-dimensional, 
compressible, nonviscous, steady flow through mixed-flow compressors 
in which the center line of the passage generates a right circular 
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cone. The two-dimensional flow pattern is considered to lie upon 
the surface of this cone. The radial-discharge compressor is a 
special case in -which the cone angle is 180°. 

The solution of two-dimensional, compressible-flow equations 
is accomplished by a numerical procedure known as the relaxation 
method. This method was first developed by Southwell (refer- 
ences 6 and 7) and has been extended to compressible-flow problems 
by Bnmons (reference 8) . It is essentially the procedures outlined 
by Emmons that are employed in the numerical solution of the dif- 
ferential equation obtained in this analysis . 

The analysis is first developed for arbitrary blade shapes 
and is later applied to logarithmic blade shapes of which the 
straight blade lying on a conic radius is a special case. A numer- 
ical example is presented with constant flow area and 20 straight 
blades operating, at a tip Mach number of 1.5. Finally, a simpli- 
fied analysis for straight blades lying on conic radii is developed 
that checks the results of the relaxation solution except for 
regions of f low near the impeller tip (and inlet) . 


AMIXSIS 
General Case 

This analysis develops a method whereby the streamlines, the 
velocity profiles, and the pressure distributions can be determined 
for steady, two-dimensional, compressible flow in centrifugal com- 
pressors with arbitrary blade shapes, varying flow areas, and fixed 
cone angles. The analysis is limited to mixed-flow compressors in 
which the center line of the passage generates a right circular 
cone with a cone angle a (fig. 1). (All symbols are defined in 
appendix A. ) The two-dimensional flow pattern is considered to lie 
upon the surface of this cone. A developed view of the conic sur- 
face is shown in figure 2. For the special case in which a is 
180°, the cone surface becomes plane and is normal to the axis of 
rotation. Such compressors (a = 180°) shall be designated radial- 
discharge centrifugal compressors. For the special case in which 
a is 0°, the cone surface becomes cylindrical and is concentric 
with the axis of rotation. Such compressors (a = 0°) are designated 
axi al-flow compressors and are not considered in this report. 

Coordinate system. - Let E and 6 be the dimensionless 
conical coordinates of a fluid particle relative to the rotating 
impeller (fig. 2). The conic-radius ratio E is defined as 
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( 1 ) 


where 

r conic radius (distance along conic element from apex of cone) 
Subscript 
T impeller tip 

The coordinate system (R, 0) rotates with the angular velocity of 
the impeller a). A particle of fluid located on this coordinate 
system is shown in figure 2. The passage-height ratio H of the 
particle in the direction normal to the conic surface (fig. 2) is 
a continuous function of the conic-radius ratio R. 


H=2L = f(R) (2) 

hr 

where h is the passage height at any conic-radius ratio R. The 
shape of the blades on the conic surface is arbitrary. 


Assumptions and limitations. - This analysis assumes that the 
flow varies only along the surface of the cone, that is, that flow 
conditions are a function of the two variables R and 0. In 
order that the assumption of two-dimensional flow on the conic sur- 
face be valid, it is necessary that the flow be uniform across the 
flow passage normal to the conic surface. In order to satisfy this 
flow condition, it is necessary that: (l) the gradient of h with 

respect to r be snail, and (2) the cone angle a (fig. 1) be 
sufficiently large. The allowable variation in a from 180° will 
depend upon the relative magnitudes of h and r and upon the 
desired accuracy. For the hypothetical limiting case in which the 
ratio h/r approaches zero everywhere alo n g the conic surface, 
the analysis is accurate for all values of a. In most practical 
centrifugal compressors, the ratio h/r is smallest near the tip 
of the impeller; therefore, the analysis is most accurate in the 
region near the tip. 

The assumption of steady flow relative to the impeller applies 
within the Impeller and for the distances upstream and downstream of 
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the impeller required to set up the boundary conditions. Stationary 
prerotation -vanes (ahead of the impeller) and stationary diffuser 
vanes (after the impeller) introduce pulsations relative to the 
rotating impeller and thus make the flow unsteady. However, these 
pulsations rapidly diminish upstream and downstream of the vanes so 
that the f low can be treated 'as steady within the region near the 
impeller provided the vanes are not too close to the impeller. 

Continuity and the stream function. - A fluid particle is 
shown in figure 2 on a developed view of the conic surface. If u 
and v are the tangential and "radial" (along the conic element) 
components, respectively, of the velocity relative to the impeller, 
then, from continuity considerations for steady flow, the following 
expression is obtained: 

i 

^ (pvHR) + ^ (puH) =0 • (3) 


where p is the weight density of fluid. 

A dimensionless stream function ^ is defined such that 


and 


=^r HR 


Po c o 




^o o 


where 

c local speed of sound 
Subscript 


(4) 


(4a) 


o absolute inlet stagnation condition 

and where the coordinate subscripts (R and 9 , in this case) refer 
to partial derivatives with respect to the coordinates. 

The continuity eguation (3) then becomes 
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Jjj ( c oPo^e^ " ( c oPo^E) “ 0 
or 

^se " = 0 

■where the double coordinate subscripts (SR and E0, in this case) 
refer to second partial derivatives with respect to the coordinates. 
Therefore, the stream function -vj/ satisfies the continuity condition. 

Irrotationality. - In the absence of viscosity, shock, non- 
uniform heat addition, and so forth, the absolute motion of the fluid 
particle is irrotational. Therefore, the absolute circulation T 
about the particle is zero and, from figure 2, 

r = 0 = jjcor^B sin ^ + u^ Rd0 

where (cor^JR sin ^ + u) is the tangential component of the absolute 
velocity. After simplification. 



Substitution of the stream function \J/ as defined by equation (4) 
gives 

1 d ( V\ 

pH / + R 3© \pHR/ 


- 'I'E ( l0 Se H )r 


- % (lose f) (6) 

E 2 \ Po J 0 


2carm sin — 


‘T 81n 2 % S 

pHE + W. 


Po c o 


which reduces to 


-^EE + X + 






dE - ^ (vdR) d0 


where the impeller tip Mach number is defined as 
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% = 


o»p sin H: 



( 7 ) 


Penalty ratio. - The general energy equation la used to deter- 
mine those terms in the differential equation (6) involving the 
density ratio p/p 0 . 

Jc P T + ^ 


(tor^ sin | + u) 2 + v 2 = JcpT 0 



( 8 ) 


where 


J mechanical equivalent of heat 

Cp specific heat at constant pressure 

T absolute temperature 

g acceleration due to gravity 

AT W temperature correction for whirl ahead of impeller 


The term 


tOnpR sin — ^(orrpR sin ~ + u^ 


S 


is the work done by the impeller upon the fluid particle if the 
fluid particle possesses no whirl (radius, r sin a/2 , times tan- 
gential component of absolute velocity) ahead of the impeller, and 
the term , a \ 

corvpRj sin -r (corpRi sin ~ + uj) 

JCpAT v sA-g L (8a) 

is the reduction in work that results from the whirl ahead of the 
impeller. The subscript I refers to impeller inlet. These same 
expressions for the work done apply in the vaneless portion of the 
diffuser, where the tangential component of the relative velocity u 
decreases (takes on larger negative values) in such a m a nn er that the 
work terms remain constant. 
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Eearrangment of equation (8) and solution for the temperature 
ratio t/t c results in 

.* -1 + Z=i 
T 0 . 2 

where 

7 ratio of specific heats 

q velocity relative to impeller 






( 9 ) 


q = 'sj u 2 + v 2 (lOj 

From equation (8a), the correction term in equation (9) for the 
whirl ahead of the impeller becomes 

-J-S = (7-1) %% [%% + ^ J (9a) 


For an isentropic process, the density ratio is related to the tem- 
perature ratio by 


1 



so that 


-*-=11 

p o 


+ 


2=i 

2 




Also, from equations 


(10) and (4), 


Po c o" 



+ 



1 

2 


( 11 ) 


( 12 ) 


Equations (11) and (12) together with the general differential 
equation (6) provide three equations with three unknowns: , 

q/o 0 , and p/p Q . The solution of these equations determines the 
steady flow of compressible fluid through oentrifugal compressors 
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with arbitrary blade shapes, with arbitrary variations in the 
pas sage -height ratio, and with a, constant cone angle. Equation (6), 
which is nonlinear, can be solved (together with equations (11) and 
(12)) by relaxation methods. 

\ 

Relaxation methods. - The relaxation solution of the differ- 
ential equation determines the value of at each point of a grid 
placed within the boundaries of the problem. Values of are 
first assumed, and the residuals Q, which result from the assumed 
values of 4^, are computed for each point of the grid by expressing 
the differential equation (for example, equation (6)) in f ini te- 
difference form with the sum of all terms equal to Q instead of 
zero. The final solution is obtained by systematically varying 
(relaxing) the interior values of 4> until the values of Q at 
each grid point approach zero. For the numerical solution of this 
problem by relaxation methods, it is convenient to transform the 
R,0 coordinates so that the arbitrary blade shapes become straight 
and parallel, thus a grid of equally spaced points between the 
blades is possible (reference 8, p. 13). 

Transformation of coordinates. - The transformation of 
coordinates is given by the analytic function 

f (z) = f(Re ie ) = !(R,e) + iq (B,e) (13) 

where z is a complex variable and the £ and q coordinates (Cartesian 
coordinates in transformed plane) are the velocity potential l in es 
(| = constant) and streamlines (q = constant) in the physical 
Re-plane for incompressible flow through the stationary impeller 
(co=0) with constant passage height (H = 1). It is convenient to 
choose the £, q coordinates to correspond to smooth (shockless) 
inlet and exit conditions with respect to the impeller blades. 

The transformation results in straight parallel lines for the 
impeller-blade boundaries as these boundaries correspond to constant 
values of the incompressible stream function q. The £, q 
coordinates are given for certain blade shapes by simple analytic 
expressions (p. 17) but, for arbitrary blade shapes, it is conven- 
ient to solve the well-known differential equations for the incom- 
pressible stream function and velocity potential by relaxation 
methods (reference 8). 

Equation (6) is now expressed in terms of the new variables £ 
and q. This change of variables is described in appendix B with 
the following result: 
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■ ~~~2 ^ T i ~ [( lo Se E )| V 1 - (• Lo Se H )q %] ( 14 ) 

where H is now a function of i and q and 

<ju velocity for incompressible flew relative to stationary 

impeller with constant blade height (HO-pl a ne) 

Ui, Vj_ tangential and radial components of velocity q-^, 
respectively 

The incompressible velocities v^, end q-j_ are related to 

the incompressible stream function q by the following equations: 


£ 

0 

1 

S3* 

(15) 

v *8 
Y 1 “ E 

(15a) 

/ 2 V* 
~lN + E s 

(15b) 


Equation (14) is the general differential equation expressed in £,q 
coordinates for the flow of compressible fluid through centrifugal 
compressors with arbitrary blade shapes, with arbitrary variations 
in the pa s sage -height ratio, and with constant cone angles. 

Also, equation (12) in terms of the transformed coordinates | 
and q (see equations (46) and (47) in appendix B) becomes 


-ea_ = i 4 2 + * 2 
p 0 °o H V e 11 

where H is a function of i and q. 


1 

2 


( 16 ) 



1015 


MCA TN No. 1744 


11 


Q 


Finite-difference equation. - In order to solve the system of 
equations (equations (11), (14), and (16)) by relaxation methods, 
equations (14) and (16) must first he changed to finite- difference 
fonn. This c h ange is accomplished with the aid of the following 
equations (reference 7, p. 19): 

, ^ 

= -k (r i - F s) 


s l< r 4- 

r u s ^ < r i + *3 - 


(17) 


where 




(f 2 + f 4 


2F) 




F any twice-differentiable function of variables £ and q 

b grid spacing 

Subscripts 

1, 2, 3, 4 four points adjacent to grid point being considered 
(F with no subscript) 

A sample grid is shown in figure 3. The grid spacing b is arbi- 
trary. However, the smaller the value of b, that is, the larger 
the number of grid points, the greater is the accuracy of the 
approximate, finite-difference equations (17). 


With the aid of equations (17), equation (14) becomes 



K 


+i ++2 + + ¥4 - <* - — 7 ~ ^ ^ 


. (ltl * ~ /w f* _ w £i\ . -A. T/,1, _ * * . 
' " p 0 y ' v l' t i ■ '5' -i 


4 y — °e p o ” *''=■© 


- W 4 - %) %] [(log© H x - lose H 3 ) v ± - (logg H 4 - log e H 2 ) Uj] 


- 5 f Hb 2 o Q . (18) 

— ^ Pr\ 

Hi 


Equation (18) Is the finite-difference equation used to compute the residual Q at each 
grid point from the estimated values of \|/' at the adjacent grid points. 


Equation (16) in finite-difference form becomes 


P<1 „ ^ 

p 0 CDtL 
O O 




(19) 


The density ratios in equation (18) oan nov be determined in the following manner: 
Prom equation (11), for fired values of BMp, curves of logg p/p 0 plotted against 
oa /o„ 0 _ can be commit ed for a c ons tant value of AT.-Al 1 -. TT-msc nurves have "hnon n "lot-bad 

1 *« 1 U W * ----- - W / ~KJ - JT 

in figure 4 for AT v /r 0 equal to zero, that is, for the case in which there is no 
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whirl in the fluid ahead of the impeller. The logarithm of 
the density ratio log e p/p 0 at each point on the grid is there- 
fore determined "by first obtaining the flow-rate ratio pq/p 0 c 0 
from equation (19) after which log e p/p 0 can he obtained from 
plots of equation (ll), such as figure 4. 

Boundary considerations. - The solution of equation (18) depends 
upon the boundary values of the stream function ty, which are 
determined from the compressor design characteristics and operating 
conditions. The boundary values of along the blade surfaces can 
be determined from the following considerations. 

The differential flow rate between adjacent streamlines is 
shown in figure 5. The radial component of the flow rate is 

pvhpTVpHRde 

and the tangential component of the flow rate is 

- puh^rvpHdR 

The differential flow rate is therefore given by 

dw «= pvhprijpEBde v puh^r^HdR 


where w is the flow rate between streamlines. From equa- 
tions (4a) and (4b), 


dw = PoCj^jjTnp (^ 0 d0 + ^ B dR) 


but 


dty = ‘vj/gdS + 


Therefore 

dw * P 0 c o h T r T 4 ^ ( 2 °) 

or 


Aw = p^hprvpAiJr 


(20a) 
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From equation (20a), tlie stream function must lie constant along 
the "blade surfaces (because the f lav rate between blades is 
constant). 


By definition, the driving face of a blade is the face in the 
direction of rotation, and the trailing face is the opposite side of 
the blade. Integration of equation (20) across the impeller passage 
from the driving face of one blade and the arbitrary setting of v 
and equal to zero along the driving face results in 

v = P 0 c 0 h T r I^ ( 21 ) 


The value of ^ at any point in the flow field is therefore a 
measure of the flow rate between the streamline at that point and 
the driving face of the impeller blade. In particular, along the 
trailing face of a blade, is constant and equal to ^ where 
the subscript t refers to the trailing face, and w is equal to 
W/B where W is the compressor (total) flow rate and B is the 
number of compressor passages (or blades). 

Equation (21) then becomes 

I = p o c o 1 ¥ r T^t ( 22 ) 

This equation can be simplified by the following considerations: 

The flow area at the impeller tip euj is given by 

a<p = BCyr^hp (23) 

where a is the included angle of the passage defined by 

O = 0 t - e d (24) 

where the subscript d refers to the driving face of a blade. 
Combining equations (22) and (23) to eliminate EhpTji results in 


¥ 

P 0 SqiCo 


*t 

0Q 1 


The boundary value of \lf along the trailing face of the blade is 
.therefore given by 
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=cpCT T 

where the flow coefficient cp is defined as 

cp = _JL_ 
<>o*T 0 o 


(25) 


(26) 


The "boundary values of in the diffuser are determined from 
the Kutta condition, which requires that the streamline along the 
"blade surface in the impeller be tangent to the blade at the tip. 
The effect of the Kutta condition on the boundary values in the 
diffuser will be discussed later in connection with the numerical 
example of this report. 

The boundary values of ^ ahead of the inlet to the impeller 
depend upon the whirl of the fluid approaching the inlet. This 
whirl determines the gradients of with respect to i and tj 
along the boundaries rather than the values of \|/ itself. If the 
flow for smooth entry is desired, this condition can be obtained 
in essentially the same way as the Kutta condition at the outlet. 
The whirl ahead of the inlet is then uniquely determined by the 
smooth inlet condition. 

Equivalent mixed-flow compressors. - As previously stated, 
this analysis applies to mixed-flow compressors with any fixed 
cone angle including the special case of the radial-discharge ' com- 
pressor in which the cone angle is 180°. However, as a result of 
the assumption in this analysis that the flow varies only along 
the conic surface, which can be developed into a plane (fig. 2), 
the solution obtained for a given fixed cone angle also applies 
to certain other mixed-flow compressors with a fewer or greater 
number of like passages with the same included passage angles, 
blade- thickness distributions, and so forth. For example, the 
included angle of the developed view of a mixed-flow compressor 
(fig. 2) is given by (assuming the blade thickness at the tip is 
zero) 


2ji sin = BQrji 

which, for a radial-discharge compressor (a = 180°) with 20 blades 
(for example), becomes 


2n = 20 G t 
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The ratio of these equations results in an expression for the mixed- 
flow cone angle a as a function of the number of impeller f low 
passages that are similar to the flow passages (20) in the radial- 
discharge Impeller 


a = 2 sin” 1 (27) 

From equation (27), the mixed-flow compressor cone angle a has 
been computed for several values of B and the results are given 
in the following table: 


Cone angle, a 
(cleg) 

Humber of passages, B 

180 (radial) 

20 

143.6 

19 

128.4 

18 

116.4 

17 


The flow in a mixed-flow compressor is therefore equivalent to 
the flow in a radial-discharge compressor with a greater number of 
blades (but with the same included passage angles a ) . Increasing 
the. number of blades in a radial-discharge compressor decreases the 
blade loading and the maximum local Mach number and increases the 
impeller slip factor. The blade loading and the maximum local Mach 
number are therefore less and the slip factor is greater for a 
mixed-flow compressor than for a radial-discharge compressor with 
the same number of blades (because 0 is less for the mixed-flow 
compressor) . 

Equivalent inward-flow turbines. - The reversal of the flow 
direction through a centrifugal compressor and the reversal of the 
direction of impeller rotation results in an inward-flow turbine. 

The general analysis presented in this report therefore also applies 
to inward-flow turbines with arbitrary blade shapes, with arbitrary 
variations in the passage-height ratio, and with constant cone 
angles. In fact, any solution obtained for a centrifugal compressor 
with smooth (shockless) entry is also a solution (with the flow 
direction and the impeller rotation direction reversed) for an 
inward-flow turbine with shockless entry. The shockless entry 
for the compressor corresponds to the Eutta condition for the tur- 
bine and vice versa. 
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Particular Case - Thin Logarithmic-Spiral Blades 

In the general case of this analysis, the differential equa- 
tions governing the flow through compressors with arbitrary blade 
shapes have been developed in terms of the transformed coordinates £ 
and q, which can be determined by relaxation methods. 

In the particular case of this analysis, a simple ana lytic 
expression that determines a particular thin-blade shape is 
considered. Such an analytic expression is desirable because it 
specifies the transformed coordinates £ and q in terms of R 
and 0 and therefore eliminates the additional relaxation solution 
otherwise required to obtain £ and q. 

Logarithmic- spiral blades. - This class of blade is given by 
the analytic function 

f(z) = (N + IK) log e z 

= (W + IE) log e (Re 10 ) 

= (N + iK) (Log e R + i0) 

= (N log Q R - K0) + i (K log e R + N0) (28) 


where N and K are constants to be determined by boundary con- 
ditions. Prom equations (13) end (28), 


£ = N(log e R - | 0) (29) 

q = N(| log e H + &) (29a) 

Any incompressible streamline (q = constant) or portion thereof 
may be replaced by an infinitely thin blade without disturbing the 
flow pattern. Equation (29a) with q constant therefore represents 
a thin-blade shape. Any number of blades may be selected corre- 
sponding to different constant values of q. 

The constants N and K are determined in appendix C from 
the boundary conditions, 
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H =*t 


when 


0 = CJq, 

and 


B = 1.0 


(from -which the zero point for 9 is at the tip of the driving 
face of the blade). Equations (29) and (29a) become 

t = cp (logg E + 9 tan 8) (30) 

t] = cp (e - tan log Q E) . (30a) 


where 8 is the blade curvature angle between conic-radius ratio E 
and tangent to blade surface (fig. 6). 


The blade shape in the B0-plane is given by equation (30a), 
which, with q constant, determines a logarithmic spiral (fig. 6). 
In appendix C, the angle 8 is shown to be constant for a given 
logarithmic spiral. Eor positive values of 8> the blades are 
curved forward (in the direction of rotation); for negative values 
of 8i the blades are curved baclcward; and for 8 equal to zero, 
the blades are straight and radial along conic radii (fig. 6) . 

Equations (30) and (30a) relate every point in the B0-plane 
to a corresponding point in the Iq-plane. Solving these equations 
for E and 0 in terms of \ and rj 

j - q tan 8 
E = e ^ 00 ® 2 P 


Q _ T| + j tan 8 
cp sec 2 8 

In terms of the particular transformed coordinates given by 
equations (30) and (30a), the differential equation for the flow of 
compressible fluid through centrifugal compressors with logarithmic- 
spiral blades becomes 


(31) 

(31a) 
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2Mj. 


-2- H 
Po 


2( j -ri tan 0) 
cp sec 2 0 

(cp sec p) 2 



+ 






- cos 0 - ^ sin 0) £( logg H)| cos 0 - (log e H)^ sin 0j (32) 

Passage-height ratio H. - Instead of the arbitrary function 
for the passage-height ratio given by equation (2), let the varia- 
tion in H be given by 


H = BP 


(33) 


where m Is an arbitrary exponent . 

If a equals Orp for all values of R, the variation in the 
flow-area ratio with R becomes, from equations (23) and (33), 


A = 


a 

®T 


BOprh 

BotpTiphjp 


= RH = tf Q+1 


(34) 


where A Is the flow -area ratio. IMs relation between A and 
R is plotted in figure 7 for several values of the exponent m. 
For m equal to zero, -the passage height 1 remains constant; for 
m equal to -1, tie passage area remains constant. 

For logarithmic -spiral blades, equation (31) is combined with 
equation (33) and the passage-height ratio H becomes 
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m(|- ti tan 3) 

H = e * se ° 2 P (33a) 

from which equation (32) reduces to 


2Mn — 

T Po 


(m +2)(£ - n tan 3) 
cp sec 2 3 


(cp sec 3)' 


% + ^ 


£-) 

H o/| 


(3*8. 


£) — Vr < V * V 

\ P °/J] cp sec 3 b 


(35) 


Straight radial blades. - For straight "blades along conic 
radii, the "blade angle 3 is zero (fig. 6) and equation (35) 
reduces to i 

(m +2)1 

^ = Hi + *nn “ % (^° Se 


- (**. £) - = t, 


(36) 


Incompressible flow. - For incompressible flow, p/p Q is 
equal to 1.0 and equation (36) becomes 


2% 


(m + 2)t 

e ? 


m 


cp 2 


+tj + - * * £ 


( 37 ) 


where Mp is now a fictions tip Mach number based upon a constant 
fictious speed of sound, which also appears in the definitions 
of |, T], and cp. 

Constant blade height. - For constant blade height, the expo- 
nent ~m is zero (equation (33)) and equation (37) reduces to 
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2i 

(38) 

Detailed solutions to this equation (in slightly altered form) have 
been obtained by a number of investigators. (For example, refer- 
ence 1.) 

Summary of equations. - Special forms of the basic differential 
equation (14) for steady flow through conical centrifugal compres- 
sors have been developed for various designs aria operating conditions. 
The equations are listed in the following table: 


Equation 

Blade shape 

Passage height 

Flow 

(32) 

Logarithmic 

spiral 

H = f (E) 

Compressible 

(35) 

Logarithmic 

spiral 

H a BP 

Compressible 

(36) 

Straight 

H = I? 1 

Compressible 

(37) 

Straight 

H » EP 

Incompressible 

(38) 

Straight 

H = constant 

Incompressible 


These differential equations express the stream function ^ 
as a function of the transformed coordinates | and rj. The equa- 
tions are changed to finite-difference form and solved by relaxation 
methods. This form is obtained with the aid of equations (17) and 
has been obtained for equation (14), the finite-difference fora of 
which is given by equation (18). The density ratio p/p Q is deter- 
mined from figure 4 or from equation (ll) with the aid of equa- 
tion (19) in which the parameter q^ is obtained from equations (10) 
and (47) or from equation (52b) in appendix C end in which the 
variable H is a given f unotion of B (whioh is related to i and 

n)« 


NUMERICAL PROCEDURE 

A detailed outline of the numerical procedures for the -relaxa- 
tion solution of compressible-flow problems is given in reference 8. 
The emphasis is placed herein on those features of the solution that 
are peculiar to the flow in centrifugal compressors. 
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The solution of each differential equation developed for the 
stream function ij / (equations (14), (32), (35), (36), (37), and 
(38)) requires special treatment depending upon the compressor 
design characteristics anfl the flow conditions that the equation 
represents. However, the numerical procedure is in many respects 
similar for all equations. The remainder of this section will 
therefore he concerned with the numerical solution of equation (36), 
which is the equation used in the numerical example of this report 
(straight thin blades along conic radii) . 

Parameters . - An inspection of the differential equation for 
the stream function (equation (36)), together with the auxiliary 
equations (7), (9a), (ll), and (25), Indicates eight parameters, 
which, together with the blade angle 6 (equal to zero in equa- 
tion (36)), specify the design characteristics and operating condi- 
tions of the compressor. 

• Design parameters: 

(1) Passage-height exponent m, which relates passage-height 
ratio H to conic-radius ratio B 

H = iP (33) 

In the more general case, H can be given as an arbitrary 
function of E. 

(2) Cone angle a, which is constant (fig. l) 

(3) Blade curvature angle 0, whioh for a logarithmic- spiral 
blade remains constant and for straight blades (along 
conic radii) is equal to zero 

(4) Included angle of impeller passage on conic surface o , 
which is measured from driving face of one blade to trails 
ing face of next blade 

(5) Inlet conic-radius ratio Bj, which is ratio of blade- 
inlet conic radius to blade-tip conic radius 

(6) Initial whirl of fluid ahead of impeller, which is des- 
ignated by ratio AT W /T 0 and whioh is determined by design 
configuration of compressor ahead of impeller and by flow 
rate (operating parameter) through compressor 
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( 7 ) 

(26) 


For a given compressor, this coefficient is easily shown to 
"be proportional to standard equivalent -flow-rate parameter 

(reference 9) 

where 

9 ratio of inlet stagnation temperature to standard sea- 
level temperature 

5 ratio of inlet stagnation pressure to standard sea- 
level pressure 

(9) Eatio of specific heats y, which determines variation 
in density ratio p/p Q as given by equation (ll) 

Boundary conditions. - The boundary values of the stream 
function ^ are determined by the parameters previously outlined 
and by the Eutta condition for tangency of flow at the blade tip. 

The various bo u ndary values of ^ are shown on the relaxation grid 
in figure 8. The manner in which these boundary values are obtained 
is summarized as follows: 

(1) The value of the stream function along the driving face of 

the blade (fig. 8) is arbitrarily set equal to zero. (See 

P. 14.) 

(2) The value of the stream function along the trailing face 
of the blade ty-j. (fig. 8) is constant and is given by 

% =cpcr T 


(25) 
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(3) The values of the stream function, along the boundaries 

ahead of the impeller-blade inlet depend (among other things) upon 
the inlet design parameters Rj and LT V /T 0 . For the particular 
case (used in the numerical examp le) of straight blades (0 = 0) 
extended to the origin, at which a point source is assumed (p. 27), 
the values of Rj and AT ff /r 0 are zero and the left bound of the 
relaxation grid (fig. 8) is located at some radius ratio R within 
the impeller passage. The value of the stream function along the 
left bound (fig. 8) is computed for this particular case from 

equation (41) . This equation has been found to be quite accurate 
for values of R = f (|) less than 0.65. (See p. 32.) 

(4) The values of the stream function along the upper bound 

of the diffuser (fig. 8) are estimated by equation (54) 

developed in appendix 3). These estimated values of are cor- 
rected by the relaxation methods to be subsequently discussed. 

The value of the stream function at the upper right comer of 
the grid 'J/j, (fig. 8) is first determined by the estimated values 
of ^ e . However, to satisfy the Eutta condition, the value of 1 l / r 
must be of such a magnitude that the flow leaves tangent to the 
blade surface at the blade tip. If the resulting relaxation solu- 
tion for the streamline configuration in the compressor does not 
satisfy the Eutta condition, a new value of \J/ r is selected, the 
estimated values of adjusted accordingly, and the solution 

repeated as many times as is necessary. 

(5) The values of the stream function along the lower bou nd 
of the diffuser (fig. 8) are, from symmetry considerations, 

'Jtfc less than ty Q . 

( 39 ) 

(6) The values of the stream function along the right bound of 
the relaxation grid are determined by and by the assumption 
that, along this bound, A^/Aq is a constant. For the numerical 
example, this assumption has been found to be accurate enough for 
values of R = f(|) somewhat greater than 1.20. 

Grid layout. - Saving determined the values of ^ at each 
point on the grid boundary, the values of at each of the 
Interior points is estimated and recorded upon the grid sheet 
(fig. 8), which should be sufficiently large to accommodate the 
calculations required by the relaxation procedure. After estimating 
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the interior values of 4*, the problem resolves itself into tvo 
parts: (l) calculation of the residual Q resulting from the 

estimated values of 4/, and (2) relaxation (elimination) of the 
residual by suitable adjustments in the estimated values of 4/. 

■Residuals, - In order to determine the residual at each point 
of the grid system, the partial differential equation (for example, 
equation (36)) is expressed in finite-difference foxm by methods 
out lin ed in the analysis. All terms of the equation are placed 
equal to the residual Q, as has been done in equation (18). The 
terms involving the density ratio p/p D are determined from fig- 
ure 4 or equation (ll) with the aid of equation (19). In order to 
determine the value of the density ratio at the boundary points in 
the impeller, it is necessary to extrapolate the values of the 
stream function 4> by graphical or numerical methods to obtain 
the velocities at the boundary. In the diffuser, however, this 
extrapolation is not required because the values of 4/ repeat 
themselves in a continuous manner from passage to passage. If the 
estimated values of 4* at all grid points are correct, the value 
of Q is zero at all points. If, however, the estimated values 
of 4> are incorrect, the values of Q are finite and may be 
positive or negative. 

Eelaxatlon. - With the residuals at each of the grid points 
computed, it remains to relax (that is, reduce) these residuals by 
suitable changes in the values of 'l'. In order to determine the 
magnitude of the required changes in 4 / , all terms of the finite- 
difference equation (equation (18), for example) are assumed to 
remain constant except the term 44/ . A change in the value of 4* 
will therefore cause a four-fold change of opposite sign in the 
value of Q. This change in the value of 4/ will also cause an 
equal change in the values- of Q at each of the adjacent grid 
points (fig. 3). These changes in 4/ and Q are recorded on the 
grid sheet as the work progresses. By continually relaxing the 
larger residuals any desired amount, the values of all residuals 
gradually approach zero. When this condition is reached, the • 
residuals are recomputed using the finite-difference equation and 
taking into account the new values of the density ratio. When 
the new values of Q have been computed, the relaxation procedure 
is repeated as often as necessary to achieve the desired accuracy. 

Accuracy. - Eo quantitative evaluation of the accuracy of 
relaxation solutions is available (reference 10, p. 176). However, 
because the computed velocities (and pressures) depend upon dif- 
ferences in the values of the stream function 4/ at adjacent grid 
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points, that Is, the small difference of large numbers, it is 
important to know the values of with sufficient accuracy to 
assure the desired accuracy for the pressure and velocity calcula- 
tions. In this report, the values of ^ were computed to the 
nearest 0.0001 compared with the maximum value of at the trail- 
ing face of the "blade of 0.1570. The resulting values of pressure 
and velocity are estimated to "be accurate within 1 or 2 percent. 

Final solution. - The streamline configuration can "be deter- . 
mined from the f ina l values of ^ at each of the grid points. 
(Streamlines are lines of constant . ) If the streamline leaving 
the tip of the impeller is tangent to the "blade at the tip , the 
estimated value of \J/ r (p. 24) is correct; otherwise a new value 
of must he selected and the entire solution repeated. 

After the correct distribution of on the grid is obtained, 
the pressure distribution can be determined from the density dis- 
tribution and the Mach number distribution can be deteimined from 
equation (19), the density ratio, and the speed of sound ratio 
c/c 0 , which is related to the density ratio by 


z=l 



From the preceding information, such quantities as the impeller 
slip (appe ndix E), the boundary- layer growth, and the blade loading 
can be estimated. 


NUMERICAL EXAMPLE 

Design and operating parameters. - A numerical example has been 
computed for the following design and operating parameters: 


Design parameters: 

(1) Constant flow area, -1 

(2) Cone angle, a, degrees 180 

(3) Straight blades along conic radii, 8 0 

(4) Included passage angle, a = C7j 2«/20 

(20 thin blades) 

(5) Inlet conic-radius ratio, Rj 0 


(6) Initial whirl of fluid approaching impeller 
inlet, AT W /T 0 


0 
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Operating parameters : 

(7) Tip Mach number, Mp 1.5 

(8) Flow coefficient, cp 0.5 

(9) Eatio of specific heats, 7 1.4 

1 

In order to insure steady-flow conditions, a vaneless diffuser 
has been selected although a vaned diffuser located far enough down- 
stream of the impeller would also be satisfactory. (See p. 5) A 
diagram of the impeller and the vaneless diffuser is shown in fig- 
ure 9. The passage-height ratio H varies in such a manner that 
the radial-flow area remains constant. The straight impeller blades 
are extended radially inward to the center where a point source has 
been assumed. A point source at the center of the impeller or a 
circular line source at some conic-radius ratio B must be assumed 
in the two-dimensional analysis in order to supply the impeller with 
the necessary flow rate. In this example, the straight blades have 
been extended to the center (Ej = 0) to obtain smooth (shockless) 
entry to the straight blades (without providing inlet whirl) . In 
practice, smooth entry to straight blades without inlet whirl is 
obtained by curved inducer vanes involving flow that is not along a 
conic surface. 

The results of the numerical example are presented in fig- 
ures 10 to 15. These figures are discussed. 

Streamlines . - The streamline configuration relative to the 
rotating impeller is shown in figure 10. The grid lines on the fig- 
ure indicate the grid spacing used to obtain the final relaxation 
solution. The streamlines are designated in such a manner that the 
value of the streamline indicates the percentage of flow through 
the passage that lies between the streamline and the driving face 
of the blade (right side of passage). For example, 20 percent of 
the flow through the passage lies between the streamline 0.2 and 
the driving face. It is interesting to note that 20 percent cf the 
flow occupies more than 50 percent of the flow area at a radius 
ratio of 0.86. The spacing of the streamlines is indicative of the 
velocities in the compressor; the smaller the spacing, the higher 
the velocity. The highest velocities therefore occur along the 
trailing face of the blade, the lowest velocities near the driving 
face. 


As a result of the irrotationality of the absolute fluid motion 
in the vaneless diffuser, the absolute tangential velocity of the 
fluid decreases as the radius ratio E increases. The tangential 
velocity of the fluid relative to the rotating E,6 coordinate 
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system therefore, decreases rapidly (takes on large negative values) 
in the diffuser and the relative streamlines in the diffuser are 
steeply sloped in the direction opposite to the impeller rotation, 
as shown in figure 10. 

Eddy formation. - For the design and operating conditions of 
this example, a wheel- type eddy forms on the driving face of the 
blade (fig. 10). This eddy is attached to the blade and rotates with 
an angular velocity equal and opposite to the rotational velocity of 
the impeller to. (The motion is not a simple rotation but a com- 
bination of rotation and deformation required to satisfy boundary 
conditions.) Such eddies form at low flow rates and result from 
the condition of irrotat iona llty relative to an absolute system 
of coordinates. 

An expanded view of the eddy is shown in figure 11. The flow 
rate within the eddy amounts to less than 2 percent of the flow 
rate through the impeller (as Indicated by the designation of the 
streamlines) . This low flow rate indicates low velocities within 
the eddy. 


It appears that, in an actual compressor under the influence 
of viscous shearing forces, the size and the position of the eddy 
might be unsteady. Udder these conditions, the pressure forces 
are not in equilibrium and the flow is unstable. It may therefore 
be desirable to eliminate the eddy by proper changes in the design 
and the operating conditions of the compressor. 

The eddy can be eliminated by the following methods: When 

the flow rate through the rotating impeller is zero, the eddy 
occupies the entire flow passage (reference 1, p. 1000). The eddy 
is therefore reduced and finally eliminated by increasing the flow 
rate through the impeller. The increase in the flow rate, however, 
is limited by choke in the compressor. In this event, it may be 
desirable to replace the eddy by a solid plug of light material 
attached to the impeller blade and disk. 

gutta condition. - Generally, several attempts are required 
to satisfy the gutta condition. The second attempt is shown in 
figure 12(a). This figure is an expanded view of the region in the 
vicinity of the blade tip. Streamline 1.0 does not leave tangent 
to the blade tip. In this case, to satisfy the gutta condition of 
tangency, it is necessary to reduce the circulation around the com- 
pressor by increasing the assumed value of ^ r . (See fig. 8 
and p. 24) 
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The third attempt to satisfy the Kutta condition is shown in 
figure 12(h). Streamline 1.0 is very nearly tangent to the blade 
tip. This solution was considered satisfactory. 

Slip. - The computed values of the slip factor (appendix E) 
as determined from the streamline configuration were 0.91 for the 
second attempt to satisfy the Kutta condition and 0.90 for the third 
(final) attempt. 

Constant Mach number lines. - Lines of constant Mach number 
relative to the impeller are shown in figure 13. A maximum local 
Mach number in the impeller of 0.64 occurs on the trailing face of 
the blade at a radius ratio of 0.68. 

The maximum relative Ifech number at the impeller tip is less 
than 0.30. The low relative Mach numbers at the tip result from 
the high tip speed of the impeller, which increases the density and 
therefore reduces the velocity. 

The somewhat irregular velocity profile at the impeller tip 
becomes quite uniform within 6 to 8 percent of the impeller radius 
beyond the tip. This rapid adjustment of the flow indicates that, 
provided the absolute leaving velocity is subsonic, vaned dif- 
fusers can be located quite close to the impeller tip without 
appreciable losses resulting from poor (unsteady) velocity distri- 
bution relative to the stationary diffuser vanes. However, in the 
presence of boundary layer and separation, the velocity profile would 
be considerably more irregular than obtained in this solution. The 
relatively close spacing of the Mach number l i nes along the trailing 
face of the blade near the tip Indicates rapidly decelerated flow, 
which is conducive to boundary- layer separation. 

Gonstant-preasure-ratlo lines. - Lines of constant pressure 
ratio are shown in figure 14. Within the Impeller, the pressures 
on the driving face of the blade are higher them the pressures on 
the trailing face. This difference in pressure across the blade 
accounts for the impeller torque. 

An expanded view of the pressure ratios in the vic ini ty of the 
blade tip is shown in figure 15. As a result of the Kutta Condi- 
tion, the pressures (and the velocities) are equal on both sides of 
the blade tip. This unloading at the blade tip decreases the 
Impeller work and therefore decreases the impeller slip factor. 
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SIMPLIFIED ANALYSIS 

The relaxation solution presented in this report is lengthy. 

It would therefore he advantageous to have a quicker, although less 
accurate, means of estimating the flow conditions within the impeller. 
In this section, a simplified analysis developed in appendix F for 
radial- di s charge and mixed-flow impellers with straight blades (along 
conic radii) is discussed. 

Telocity distribution. - The simplified analysis is based on 
the assumption that the tangential component of the velocity rela- 
tive to the impeller is zero at all radii within the impeller. This 
simplified analysis includes blades of varying thickness provided 
the variation is not sufficiently great to result in appreciable 
values of the tangential component of velocity to the impeller u. 

The irrotationality equation (5) reduces to 

a 1 civ 

sto 5 * S Si 

which, when integrated between limits, and so forth, becomes 

♦ 

= -T + 2BMp 0 (40) 

o o 

where v^ is the velocity along the driving face of the blade 
(0 =0), which is determined from continuity considerations in the 
following paragraph. 

Streamline distribution. - From continuity considerations and 
from the velocity distribution given by equation (40), the follow- 
ing expression is obtained for the streamline distribution 
(appendix F)i 
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where the constant term AT w /r 0 is evaluated by equation (9a). The 
velocity ratio v^/c D is obtained "by equation (41) from the condi- 
tion that 


ii = 'i't = <P0T 

when 


0 = 0 

When the velocity ratio v^/c D is determined, equation (41) 
expresses the streamline distribution as a function of the coordi- 
nates R,0 and the tip' Mach number Mp. 

Discussion. - The velocity ratio along the driving face of the 
blade v^/c Q tea been computed from equation (41) for the same con- 
ditions used in the relaxation solution of this report. The results 
of this computation are shown in figure 16 and are compared with the 
results obtained from the relaxation solution. The agreement 
between the relaxation solution and the approximate solution is 
satisfactory up to a radius ratio of about 0.85. For radius ratios 
greater than 0.85, the assumption that the u component of the 
velocity, and its derivative, may be neglected is no longer valid 
because appreciable changes in the u component near the blade tip 
are required, from momentum considerations, to unload the blade. 

In comparing the flow characteristics across the passage 
obtained from the relaxation solution and the approximate solution, 
the radius ratios 0.675 and 0.855 shown in figure 16 will be used. 

The eddy formed on the driving face of the blade (fig. 10) has 
two stagnation points on the blade. The velocity along the driving 
face of the blade is zero at radius ratios corresponding to these 
stagnation points (fig. 13). The smaller radius ratio can be pre- 
dicted by the approximate solution from the plot of v^/c o in fig- 
ure 16. The estimated value of the radius ratio is 0.703 compared 
with a value of 0.713 obtained by the relaxation solution. It 
therefore appears that the presence of the eddy can be predicted 
with fair accuracy provided the eddy occurs at a radius ratio for 
which the simplified analysis is valid, that is, for a radius ratio 
at which the u component of the velocity may be neglected. 

The velocity distribution across the passage has been computed 
from equation (40) with the computed values of v^/c Q given in 
figure 16. This velocity distribution is plotted in figure 17 for 
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radius ratios of 0.675 and 0.855 and is compared with the corre - 
sponding velocity distribution obtained by the re laxa tion solution. 
The velocity distribution at a radius ratio of 0.675 is nearly the 
same for both solutions. At a radius ratio of 0.855, the agreement 
is still satisfactory although the slopes of the two curves have 
begun to deviate. This deviation between the solutions results 
from the blade unloading. 

The stre amline distribution across the passage has been com- 
puted from equation (41) with the computed values of v^/c Q given 
in figure 16. This stream line distribution is plotted in figure 18 
for radius ratios of 0.675 and, 0.855 and is compared with the corre- 
sponding stre amlin e distribution obtained by the relaxation solution. 
The fact that the streamline distribution at a radius ratio of 0.675 
is nearly the same for both solutions indicates that the boundary 
values of 'l' at the left bound (fig. 8) can be computed by 

the approximate method for radius ratios at least as high as approx- 
imately 0.65. At a radius ratio of 0.855, the streamline distribu- 
tion, although satisfactory in a qualitative manner, has begun to 
deviate appreciably from the more rigorous relaxation solution. 


Pressure ratio. - The pressure ratio in the impeller passage 
is given by (appendix F) 



where v^/cq is determined from continuity considerations as 
described on page 31. 


(42) 


The pressure ratio across the passage has been computed from 
equation (42) with the approximate values of v i /c Q given in fig- 
ure 16. This pressure distribution is plotted in figure 19 for 
ra dius ratios of 0.675 and 0.855 and is compared with the corre- 
sponding pressure distribution obtained by the re lax ation solution. 
The pressure distribution at a radius ratio of 0.675 is nearly the 
same for both soluti ons and at a radius ratio of 0.855 the agree- 
ment is satisfactory. 

Blade loading, - The blade loading at any given conic-radius 
ratio is given by the difference in pressure between the driving 
face and the trailing face of the blade. From appendix F, this 
difference is given by 
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Pd ~Pt 

“p7~ = ~Sr 

■“o 


(43) 


The blade loading has been computed from equation (43) and is 
plotted in figure 20 together with the corresponding curve obtained 
from the relaxation solution. The two curves are in good agreement 
for radius ratios less than 0.825. At higher radius ratios, how- 
ever, the blade unloads for the relaxation solution whereas the 
blade loading continues to increase for the approximate solution. 

In order to correct equation (43) for the blade unloading, the fol- 
flowing relation is assumed: 


p a -p t 27 EHt„_ „ ^ , 



where the exponent x is determined by the slip factor p from 
considerations given in appendix F 


x 


= 

1- p. 


(45) 


In the numerical example, the value of p obtained from the 
relaxation solution is 0.90, for which, from equation (45), the 
value of x is 18.0. Equation (44) for the blade loading there- 
fore becomes 

?d-Pt 
Po 


27EHp 


cpa T (l-R 18 ) 


HR 


(44a) 


CONCLUSIONS AND SUMMARY OF RESULTS 

A general method of analysis has been developed for two- 
dimensional, steady, compressible flow in centrifugal compressors 
with arbitrary blade shapes, arbitrary variations in the passage 
height, and with fixed cone angles (right circular cones generated 
by the center line of the flow passage). From this analysis the 
following conclusions can be drawn: 
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The solution obtained for a given cone angle also applies to 
other cone angles (that is, other mixed-flow compressors) with a 
fewer or a greater number of passages but with the same included 
passage angle, and so forth. Mixed -flow compressors with the 
same number of flow passages as radial -discharge compressors and 
therefore with, smaller included passage angles have lower peak 
blade loadings and lower maximum relative velocities than the 
corresponding radial -discharge compressors. 

The general analysis also applies to inward-flow turbines, 
lhe solution obtained for a centrifugal compressor with smooth 
(shockless) entry is also the solution (with the flow direction and 
rotation reversed) for an inward -flow turnbine with the same design 
characteristics (that is, the same rotor) and with, smooth (shockless) 
entry. 

kram a numerical example, of a radial-discharge compressor with 
20 straight blades lying on conic radii, constant flow area, an 
i n let radius ratio of zero, a tip Mach number of 1.5, and a flow 
coefficient of 0.5, the following results were obtained: 

1. A wheel -type eddy formed on the driving face of the impeller. 

2. Hie velocities and pressures at the impeller tip were 
reasonably uniform and any nonuniformity in the flow leaving the 
impeller rapidly adjusted itself. 

3. Hie maximum, local Jfach number relative to the impeller was 
0.64 and occurred along the trailing face of the blade at 68 per- 
cent of the tip radius. 

4. Hie computed slip factor was 0.90. 
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In addition, a simplified analysis for straight blades lying 
on conic radii is presented which can be used to determine the 
streamlines, pressure distribution, and velocity profiles within 
the impeller except near the tip (and inlet). 


Lewis Flight Propulsion Laboratory, 

national Advisory Committee for Aeronautics, 
Cleveland, Ohio, August 18, 1948. 
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APEEHDIZ A 
SYMBOLS 

The following symbols are used in. the analysis: ' 

A flow-area ratio, a/a^ 

a flow area, normal to conic surface 

B number of blades (or passages) 

b grid spacing (fig. 3) 

c local speed of sound 

Cp specific heat at constant pressure 

F any tw ice-differentiable function of two variables 

g acceleration due to gravity 

H passage-height ratio, h/hrji 

h passage height, normal to conic surface 

J mechanical equivalent of heat 

K constant 

( 0 r„ sin § 

Mj impeller tip Mach number, — ± — 

c o 

m passage-height exponent 

N constant 

p static pressure 

Q residual 

q velocity relative to impeller, + v^ 

R conic-radius ratio, r/r T 

r conic radius (distance along conic element from apex of cone) 

T absolute temperature 
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u tangential component of velocity relative to impeller (positive 
in direction of rotation) 

v radial (along conic element) component of velocity 
V compressor flow rate (total) 
w flow rate between streamlines 
i exponent for blade-loading correction 
z complex variable 

a cone angle (fig. l) 

P blade curvature angle (fig. 6) 
r absolute circulation 

7 ratio of specific heats 

A finite increment 

8 ratio of inlet stagnation pressure to standard sea-level 

pressure 

T] transformed coordinate ( incompressible stream function) 

0 angle, radians (fig. 2) 

9 ratio of inlet stagnation temperature to standard sea-level 

temperature 

M- ' slip factor 

£ transformed coordinate (incompressible velocity potential) 

p weight density of fluid 

a included passage, angle, (9^. - 0 d ) 

T impeller torque 

cp flow coefficient, 

’ Po a T°o 

\J / compressible stream function 

u) impeller angular velocity 



38 


MCA Err Ho. 1744 


Subscripts: 

av 

d 

e 

f 

I 

1 

Z 

o 

r 


average 

driving face of blade (fig. 8) 

upper bound of diffuser (fig. 8) 

lower bound of diffuser (fig. 8) 

impeller inlet 

incompressible flow 

left bound (fig. 8) 

absolute inlet stagnation condition 

right bound (fig. 8) 


S 

T- 

t 

w 

R, 0, | , and 


ER, 00, R0, ee, 
ipi, and 

1, 2, 3, and 4 


simplified solution 
impeller tip 

trailing face of blade (fig. 8) 
whirl ahead of impeller 

partial derivatives with respect to R, 0, | , arid 
t] } respectively 

second partials with respect to R, 0, R and 0, 

T), and | and tj, respectively 

grid points adjacent to point in question (fig. 3) 
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APPENDIX B 


TRANSFORMATION OF GENERAL DIFFERENTIAL EQUATION 
FOR FLOW IN CENTRIFUGAL COMPRESSOR 
FROM R,0 TO |,ij COORDINATES 

The conformal transformation (except for singular points) freon 
the R0- to the Iq-plane, where' | and T| are the Telocity potential 
and stream function, respectively, for incompress ihle flow through 
the stationary impeller (to = 0) with constant passage height (H «= l), 
is given "by the analytic function 


f(z) o ffRe 1 *) o £ (R,0) + iri (R,0) (13) 

where £ and tj are functions of R and 0. Therefore, 


f r 

= F^ | s + 

F H 

^R 


f rr 

- F t e fe 2 

+ ; 


1 >s + 

F 0 

- F e £ e + 




F 00 


+ : 

2F* 

h 

h le + 

velocity potential 

S 

and the 


> (46) 


'00 


hy the Cauchy-Riemann differential equations 


i b R R 


Ui = T 


" % 


(47) 


After equations (46) and (47) are combined, the first three 
terms of the right side of equation (6) become ' 
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^B ^00 

B B 2 





% ri ee\ 

b b 2 ; 


But, "because f(z) is analytio, 


’lfi ^00 
^BR + — + ~2 
B B^ 


and from equation (47), 



where q^ is the velocity for incompressible flow through the station- 
ary impeller ((o*= 0) with constant blade height (H « 1). Therefore, 


'I'p 

*BH+— + 
B 





(48) 


The fourth term of the right side of equation (6) becomes 

Uog e H) fi «= d°Se H >| If + < lo 8e E \ %j 

* -*n ®l) (l0g e E) t V i ' (1 ° S « H)l 1 ^ 


(49) 


where H is now a function of | and T| instead of B alone. 


The last two terms of equation (6) in like manne r become 




_£. \ H -s- + <*2 


T^VV H L V 4 

_ „ Btl<M (48) to (SO) , (6) 1>e0aM8 

p^fij&r ocuibinii 1 ® *L 

*5 A , .* w ■ n (^ *V *" ^ ^ 


4i® <>o 


V 




) % 


(50) 


y 

£ 


(14) 


*>■ 

H 


♦ o£ SI TO^EI 
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APPENDIX C 

DIFFERENTIAL EQUATION FOR FLOW THROUGH CENTRIFUGAL 
COMPRESSORS WITH LOGARITHMIC-SPIRAL BLADES 
Consider the analytic function. 


f(z) «= (N + IK) log e z 

« (N + IK) log 0 (Re 1 ®) 

» (N + IK) (log e B + i0) 

* (N log 0 R - K0) + i(K log Q R + N0) (28) 

where N anri K are constants to "be determined hy the boundary con- 
ditions. From equations (13) and (28), 

| * N(log e R-|e) (29) 

and 

q «= N(| logg B + 0) (29a) 

The blade shape in the RS-plane is given by equation (29a), 
which, with q constant, determines a logarithmic spiral (fig. 6). 

Tn order to determine the constants N and E, equation (29a) 
is differentiated with q constant. 


0 


E dR 
N R 


“ + d0 


K Rd0 
H = “ dR 


But, from figure 6 
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P 


tan* 1 SW 
dB 


viier© p is the angle "between "blade surface and conic-radius ratio B 
on Be-plane (and on conic surface). Therefore, 

“tr 

== s - tan 8 « constant 
N 

and equation (29a) "becomes 

t] « N(0 - tan p logg B) (51) 

In order to determine the constant N, let 

* - 

i 

when 


8 * (Jj 

and 


B » 1.0 

(from vhioh the zero point for 8 is at the tip of the driving face 
of the blade). Therefore,' from equation (51), 



which, from equation (25), is equal to the flow coefficient cp and 
the final expressions for £ and tj become 


£ = cp ( log 0 B + e tan p) (30) 

T| = cp(0 - tan p log e B) (30a) 

These equations relate points in the Bg-plane to corresponding points 
in the | T| -plane. Solving equations (30) and (30a) for B and 8 in 
terms of £ and tj gives 
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R *= e 


l Z -H .to n . P 

2 

cp sec p 


0 = 


_ ti + l ten ft 


cp sec ft 


From equations (15), (30a), and (31), 

t| tan ft - j - 
* cp tan ft e ^ 090 ^ 


and 


so that 


v.^ *= cpe 


n tan ft-£ 

2 a 
cp sec ft 


*i “ V u i 2 + v l 2 


(31) 

(31a) 


(52) 


(52a) 


( 10 ) 


-cpsec ft e 


* to*. &r i 

2 0 

cp sec ft 


(52b) 


From equations (14) and (52) to (52b), the differential equation for 
the flow of compressible fluid through centrifugal compressors with 
logarithmic-spiral blades becomes, in terms of the transformed 
coordinates | and tj, 

2(g- t) tan ft) 

2 

cp sec 3 / \ / 

t * . : — ■ ' + '"w - ♦eH* th ' l l08e t 


( cpsec ft) 2 


- (ty cos ft sin ft) £( log 0 H)| cos ft - (log e H)^ sin ft 


] 


ojt\ 


(32) 
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APPENDIX D 


PROCEDURE FOB ESTIMATING BOUNDARY VALUES OF STREAM 
FUNCTION IN VANELESS DIFFUSERS 

Estimated ‘boundary values of the stream function in vaneless 
diffusers can he obtained by assuming, as a first approximation, 
that the flow in vaneless diffusers is one dimensional, that is, a 
function of R only. From the conservation of the moment of momentum. 


r sin 




r T sin | sin 


where n, the slip factor, is defined in appendix E. Therefore, 
solving for u, dividing by c D , and expressing in terms of R 



(53) 


In general, the slip factor p depends upon the blade shape ft-nS 
upon the operating conditions of the compressor. Combining equa- 
tions (4a) and (53) gives 



(54) 


This equation gives the variation in the stream function with the 
radius ratio R for a constant value of 9. Because the density 
ratio is also a function of R, equation (54) is solved by numerical 
point-by-point methods. In order to obtain this solution, it Is 
first necessary to Enow the density variation with R. 

The density ratio is given by 




(11) 


v 
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where 


©‘ • (t)‘ • ft) 1 

which, from equations (26) and (53) and from continuity consider- 
ations, assuming a^ = 2« sin hj r^, 



The upper or lower boundary values of the diffuser \J/ 0 or ^ 

(fig. 8) are therefore estimated from equations (54) and (55) using 
assumed values of the slip factor. These estimated boundary values 
are not generally the oorreot values and provision must be made in 
the relaxation solution for their correction in the same manner 
that any other estimated values of are corrected. 
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APPEHDIX E 


PROCEDURE FOB COMPUTHKJ IMPELLER SUP FACTOR 

The Impeller slip factor is defined as the ratio of the average 
absolute tangential velocity of the air as it leaves the impeller 
tip to the tip speed of the impeller. 


slip factor o n * 


(tor T sin f -h u) ay 


cor T sin S 


1 + 


^av 

“ r T <*> f 




«T 


(56) 


The average value of the tangential velocity relative to the impeller 
at the impeller tip is given by 




P6i 


9, 


pvh^ r T d0 


Pq c o *T r T 

- I 

B o 



d0 


which, from equations (22) and (25), reduces to 




1 

5 


( 57 ) 
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Equation (58) gives a weighted average value of u/c 0 . This 
weighted average value of u/c c is also equal to the unweighted 

P e t 

average value cf u/cq, which is equal to — I /Ji-A d8. f act 

°T \°oJ 

j0 d 

can he shown from considerations of the conservation of angular 
momentum in the diffuser, which is based upon the weighted average 
value of u/c 0 and from considerations of constant absolute cir- 
culation in the diffuser, which is based upon the unweighted* aver- 
age value of u/c c . Combining equations (56) and (57) results in 
the following expression for the slip factor: 


H = 1 + 


ol - 

'PMr UoApoJUJ V 

Jo 


6 - dr 


(58) 


The value of the integral is obtained from the area under a curve 
of the integrand evaluated at the impeller tip as a function of the 
9 - 8 d 

angle ratio (fig. 2l).- The value of the integrand is 


obtained from the relaxation solution. 
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APPENDIX F 


I 

SIMPLIFIED ANALYSIS FOE EADIAL-DISCEAR® AND MIZED-FLOW 
CENTRIFUGAL COMPRESSORS WITH STRAIGHT BLADES 


ALONG CONIC RADII 


Velocity distribution. - This simplified analysis is "based on the 
assumption that for centrifugal compressors with straight blades along 
conic radii the tangential component of the velocity relative to the 
impeller is zero at all radii within the impeller. The irrotationality 
equation (5) therefore reduces to 


2oxrvp sin ^ 



dv 


which, when integrated with respect to 6 between limits, becomes 

v = v^ + Stttnjj sin ^ BO 

which, in turn, divided by the speed of sound at inlet stagnation 
conditions becomes 



+ 2RMp 6 


(40) 


where , v^ is determined from continuity considerations in the fol- 
lowing paragraph. 

Streamline distribution, - Neglecting the u component of 
velocity, the continuity equation becomes (fig. 2): 


dw = pvhrjTiVjiHRd0 


The density p is obtained from equation (11) with q/o Q equal to 
v/o G and the velocity ratio v/c 0 is obtained from equation (40), 
so that 
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dw 


VoVi 


= HR <1 


+ 2zl 


(BMp)' 




_ _JZV ( — + 2RMp0) de (59) 


•where the constant term AT^/Tq is evaluated by equation (9a) . 
Integrating the right side of equation (59) between the limits 0 
mifi q nnd noting from equation (20) that the left side of equa- 
tion (59) is equal to d 1 ^, -which is integrated between 0 and \|/, 
gives 



Equation (41) expresses the streamline distribution as a function 
of the coordinates E,0 apd the tip Mach number Mji. The velocity 
ratio v d /c Q is obtained by equation (41) from the condition that 


when 


^ = ^ a T 


9-0 

Pressure ratio. - The pressure ratio is directly determined 
from the temperature ratio by the thermodynamic relation 
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The temperature ratio is given by equation (9) in -which the velocity 
ratio q/c Q is equal to v/c 0 given by equation (40) so that 



where v^/o Q is determined from cont inui ty considerations 
viously described. 


(42) 


as pre- 


Substituting equation (42) into equation (41) results in 


* 


HR / P d _£_\ 

2 yB 2 Aji ^p Q ’ p J 


(60) 


The stream function and the pressure ratio p/p Q are there- 
fore related by the simple equation (60). 


‘ Blade loading. - The blade loading at any given radius ratio 
is given by the difference in pressure between the driving nnfl 
trailing faces of the blade. From equation (60) this difference 
is 


P d - P t 
~o 


+t 


which, from equation (25), becomes 



27EHn 

IT *°* 


(43) 


The blade loading has been computed from equation (43) and is 
plotted against the conic-radius ratio B in figure 20 together 
with the corresponding curve obtained from the re laxa tion solution. 
The two curves are in good agreement for radius ratios less than 
0.825. At higher radius ratios, however, the blade unloads for the 
relaxation solution, whereas the blade loading continues to increase 
for the approximate solution. Xn order to correct equation (43) 
for the blade unloading, the following relation is assumed, 
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p d -p t 27HH,, 


HR 


cpa g, (1- R 1 ) 


(44) 


•where the exponent x is determined "by the slip factor from the 
following considerations: 


The impeller torque t is obtained from 
equation 


dT = Bp 


' p d“ 


CL 

hr sin - 
2 


which expressed in terms of the dimensionless 
R and combined with equation (44) results in 


the differential 
dr 

coordinates H and 


dT = Bp 0 2yMj7 cpOjh»piVp 2 sin ~ Rtl-R 1 ) dR 

Integrating between the limits of R equal to 0, assuming the 
blades extend to the origin, and R equal to 1.0 gives 


T = Bp o 7^0 tVt 2 Sin | (xfg) (61) 

The torque can also be obtained from an expression for the 
impeller power, which includes the slip factor p. 


cot = ji ^cortp sin --Q 2 W 

which, after having been combined with equations (7), (23), and 
(26), the terms being rearranged, becomes 

t = Bp 0 7 ^cpa T h T n r 2 sin | (p) (62) 

Prom equations (61) and (62), the exponent x is obtained as e 
function of the slip factor p. 
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The value of |i obtained from the relaxation solution is 0.90 
r vhioh, from equation (45), the value of x is 18.0. Eq.ua- 
ion (44) for the blade loading therefore becomes 





(44a) 
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Figure 2. - Fluid particle on developed view of conic surface. R, e, and H, dimensionless coordinates 
relative to Impei ler; u and v, tanoentlal and radial components of velocity relative to Impeller, 
respectively. 
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Figure 3. - Sample grid showing grid spacing b and numerical 
subscript convention for adjacent grid points. 







Natural logarithm of density ratio, log e p/p 
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Figure 4. - Natural logarithm of density ratio log e p/p 0 as function of flow-rate ratio 
pq/p 0 e Q for various values of RMy. Equation (Jl), ratio of specific heats equal to 
1.4; whirl ahead of impeller equal to zero. 
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Figure 7. - Flow-area-ratio variation with conic-radius 
ratio for several values of exponent m. Blade thick— 
aess assumed to be a constant percentage of passage 
width at all conic-radius ratios. 
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Figure B. - Relaxation grid In transformed coordinates showing boundary values 

of v discussed In text. 
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Figure iO. - Relative streamlines for compressible flow through centrifugal compressor. 
Streamline designation indicates percentage of flow through passage between streamline 
and right side of passage. Passage angle, 18 degrees; Impeller tip Mach number Mf, 1.5; 
flow coefficient <p, 0.5; constant flow area (m, -1.0). 
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Flours II. - Vortex on driving face of blade. Relative streaallnes 
for every 0.2 percent of flow throuoh lapel ler passage. Passage 
angle, 18 degrees; lopeller tip Uaeh nuaber Uj, 1.6; flow 
coefficient ?, 0.6; constant flow area (a, -1.0). 
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G . ; 


Relative streamlines, 



H f 1 1 1 1 h h 1 — 

2.0 1.5 1.0 .5 0 -.5 -1.0 -1.5 -2.0 

Angle, deg 

(a) Second attempt to satisfy Kutta condition of flow tangency at blade tip. 


Figure 12. - Expanded view of relative streamlines in region of blade tip. Passage 
angle, 18 degrees; impeller tip Mach number My, 1.5; flow coefficient y, 0.5; 
constant flow area (m, -1.0). ' 
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Relative streamlines, y/f t 



Angle, deg 

(b} Third attempt to satisfy Kutta condition of flow tangency at blade tip. 

Figure 12. - Concluded. Expanded view of relative streamlines in region of blade 
tip. Passage angle, 18 degrees; impeller tip Mach number My, 1.5; flow coefficient 
<?, 0.5; constant flow area (o, -1.0). 
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CO 


Relative 
Mach number 



Figure 13. - Lines of constant Mach number relative to impeller. Passage angle, 
18 degrees; impeller tip Mach number My, 1.5; flow coefficient ?, 0.5; constant 
flow area (m, -1.0). 
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Pressure ratio, p/p 0 



2.0 1.5 1.0 .5 0 -.5 -1.0 -1.5 -2.0 

• Angle, deg 


Figure (5. - Expanded view of pressure-ratio lines in vicinity of blade tip. Passage 
angle, 18 degrees; Impeller tip Mach number My, 1.5; flow coefficient ?, 0.5; 
constant flow area (m, -1.0). 





Velocity ratio at driving face of blade 
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.65 .70 .75 .80 .85 .90 .95 1.00 


Conic-radius ratio, R 

Figure 16. - Comparison of velocity ratio along driving face 
of blade for approximate and relaxation solutions. 
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a 


Figure 17. - Comparison of velocity distribution 
across passage for approximate and refaxation 
solutions. 




St ream function 


NACA TN No. 1744 







Pressure ratio, p/p 
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Figure 19. - Comparison of p res s u r e- rat i o distribution 

across passage for approximate and relaxation 
sol ut ions. 



Blade 1 oad Ing , (p d -p t ) /p 
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Conic-radius ratio, R 


Figure 20. - Comparison of blade loading for approximate and relaxation 

sol ut ions. 
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Angle ratio across passage, 0-0 d 

0 T 

Figure 21. - Variation in integrand (for computing 
impeller slip) across passage at impeller tip. 
Second and third attempts to satisfy the Kutta 
condition. 



